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On Calculation of Vector Spin Chirality for Zigzag Spin Chains 

Kouichi Okunishi 
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We calculate the vector spin chirality for S — 1/2 zigzag spin chains having U(l) symmetry, 
using the density matrix renormalization group combined with unitary transformation. We then 
demonstrate the occurrence of the chiral order for the zigzag XY chain and discuss the asso- 
ciated phase transition. The results are consistent with the analysis based on the bosonization 
and the long distance behaviour of the chirality correlation function. For the S — 1/2 zigzag 
Heisenberg chain in a magnetic field, we also verify the chiral order that is predicted by the 
effective field theory and the chirality correlation function, and then determine its magnetic 
phase diagram. 
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1. introduction 

Recently the vector spin chirality in low dimensional 
quantum systems has been attracting much attention. 
In particular, the 5 = 1/2 zigzag chains have been pro- 
viding considerably interesting physics associated with 
the ordering of vector spin chirality at the zero tempera- 
ture; 1-4 ) A strong quantum fluctuation breaks the classi- 
cal helical order and the low-energy physics is described 
by disordered ground states such as the Tomonaga- 
Luttinger(TL) liquid state and dimerized state. However, 
a "hidden" vector chiral order can be realized in the anti- 
symmetric sector, accompanying the spontaneous break- 
ing of Zi symmetry originating from the double-chain 
nature of the zigzag lattice. In Refs.f 5-7 )], it was also 
shown that the vector chiral order is also induced for 
the zigzag Heisenberg spin chain by a uniform magnetic 
field. An interesting point in the chiral order of the zigzag 
chains is that it can play an important role in the exper- 
imental situation of multiferroics through the spin-orbit 
coupling, 8 ) although the chirality itself may not be de- 
tected by uniform magnetization measurement. 

The vector spin chirality we discuss here is defined 

as 

kij — [Si x 5j] z , (1) 

on the bond of the i- and j-th sites, where S is the 
5=1/2 spin operators and the spin quantization axis 
associated with the U(l)-rotational symmetry is in the 
z-direction. 9 ' 10 ) This quantity is also mentioned as a spin 
current on emphasizing the quantum nature of the U(l)- 
symmetric spin chain: jy = J kij = —iJ/2{SfSj — 
S~Sj~), where J is the exchange coupling. A theoreti- 
cal difficulty in investigating the chiral order, particu- 
larly in numerical approaches, is that, in the convention 
of the Pauli matrix, all matrix elements of eq. (1) are 
pure imaginary, while the wavefunction can be repre- 
sented by real numbers, since the spin chain Hamiltonian 



n+1 




n n+2 J 2 

Fig. 1. Zigzag lattice; n=even and odd sites form the double 
chain structure. 



is a real symmetric matrix. This implies that an exter- 
nal field conjugate to the chirality is exactly zero even 
at the level of numerics, and then the expectation value 
of eq. (1) must be zero even in the occurrence of a spon- 
taneous order of the chirality. In the previous numerical 
analysis, accordingly, the properties of the spin chirality 
were extracted from the long-distance behavior of the 
chirality-chirality correlation function. 3, 6 ' 7 ' n ' However, 
the groundstate expectation value of the chirality opera- 
tor itself is clearly desired for the precise analysis of the 
chiral order. 

In analogy with the Josephson current in supercon- 
ductivity, a phase of spin, which is associated with a spin 
rotation angle, can be expected to be crucial for spin 
current. This suggests that a certain spin rotation makes 
the spin chirality more clearly observable and enables us 
to directly calculate the spin chirality. In this paper, we 
actually introduce a unitary transformation that makes 
the Zi symmetry associated with the spin chirality more 
visibly, and then calculate the spin chirality order pa- 
rameter, using the density the matrix renormalization 
group (DMRG). 12 -* On this basis, we discuss the various 
properties associated with the chiral order for zigzag spin 
chains. 

The 5=1/2 zigzag chain is one of the most essential 
ID quantum spin systems having the frustration effect; 
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its Hamiltonian is given by 

H = Hi+H 2 (2) 

with 

n x = j\ $>*s* +1 + sy n si +l + as*s* +1 ], (3) 

n 

h 2 = j 2 ^Ks: +2 + SlS y n+2 + AS*S* +2 ], (4) 

n 

where Ji and J 2 denote the nearest and next-nearest 
neighbor couplings, respectively, and A is the anisotropy 
parameter. As in Fig. 1, the n =even and odd sites form 
the double-chain structure without losing the transla- 
tional invariance. In this paper, we basically consider 
Ji > Ji > and thus introduce the notation a = J1/J2 
for simplicity. The symmetry of the zigzag spin chain 
and vector chiral order was analyzed in Ref.[ 4 )]. Here, 
we would like to mention that the Z 2 symmetry corre- 
sponding to the parity between n=even and odd chains 
is important for the vector chiral order. 

This paper is organized as follows. In §2, we intro- 
duce the unitary transformation corresponding to the 
7r/2 rotation of the spin phase, which makes the Z 2 sym- 
metry of the system more transparent. As a similar but 
different approach, we also introduce phase modulation 
into the system, which can be regarded as a field conju- 
gate to the chirality. In §3, we calculate the expectation 
value of the chirality operator for the zigzag XY chain 
(A = 0) in the bulk limit, using the infinite system size 
method of DMRG accelerated by the recursion relation 
for the wavefunction. 13 ) We then find that the chirality 
order parameters basically agree with the previous re- 
sults based on the bosonization and the analysis of the 
correlation functions. We also discuss the chiral order in 
the context of phase modulation in §2. In §4, we calculate 
the spatial distribution of the local chirality for the zigzag 
Heisenberg chain in a magnetic field. We then show that 
the chiral order is actually induced for a certain region 
of the magnetic field and a, which is also consistent with 
the previous weak coupling theory 5 ) and the correlation 
function analysis. 6 ^ We then find that the transition be- 
tween the chiral-ordered phase and the two-component 
TL liquid 14 ^ may still appear, and finally determine the 
magnetic phase diagram in addition to the phase bound- 
ary obtained previously 15 ' 16 ) . In §5, we summarize the 
results. 

2. chirality and spin phase 

As was described in the introduction, the main dif- 
ficulty in the numerical calculation of the chirality order 
parameter is in that the Hamiltonian (2) is a real sym- 
metric matrix, which is not suitable for dealing with the 
phase of spins; the spin current is caused by the phase 
difference between n =even and odd sites, which is as- 
sociated with the Z 2 symmetry breaking of the sign (di- 
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rection) of the current. We thus try a spin rotation to 
control the relative phase of the spins on even and odd 
sites. 

2.1 unitary transformation 

Let us write the spin rotation operator around the 
z-axis at the site n as Rn(0) = e l9S ™, so that R„(9) : 
S*+ — ► e* e S+ and S* — > 5*. Then we define a unitary 
operator as, 

U Z (^)= R nfy (5) 

nGcvcn 

which rotates spins at n =cvcn sites by 7r/2 around the 
z axis. Then Hi is transformed into 

(— )" 

Hi = ^ —^—(Sn^n+l — &n &n+l) + ^^n^n+l) (6) 
n 

while H 2 is unchanged. Note that, if 9 = —ir/2 rotation 
is employed instead of eq. (5), the sign at the XY term 
in Hi is inverted. Thus, the sign of the spin rotation is 
closely related to the Z 2 symmetry of the spin chirality. 

The original Hamiltonian (2) is invariant under the 
one-site translation operation: n — > n + 1, but a naive 
translation of the transformed Hamiltonian (6) changes 
the sign in front of the XY term, which implies that the 
explicit translational invariance is hidden by the unitary 
transformation. However, the translation combined with 
an exchange of the label of even and odd sites preserves 
cq. (6) including the sign of the XY term. In the repre- 
sentation of eq. (6), the other symmetry operations pre- 
sented in Ref. [ 4 >] are also preserved. Then an essential 
point is that the Z 2 symmetry in the Hamiltonian be- 
comes more explicit, which is expected to be suitable for 
numerical treatment. 

According to (5), the chirality operator is also trans- 
formed into 

(-)" 

«n,n+l = ^— {Sn S n +1 + &n ^n+l)- CO 

As can be seen for eq. (6), the local chirality operator (7) 
is invariant for the one-site translation combined with 
the exchange of the label of even and odd sites. Then an 
important point is that the sign of eq. (7) is assigned 
staggeredly, so that the sign of the uniform operator 
Sn^n+i + Sn^n+i snou ld be broken staggeredly to pro- 
duce a uniform chiral order. In addition to the above, we 
just note that the form of the chirality operator and the 
XY term in Hi exchanges with each other in this unitary 
transformation, which may interestingly suggest a kind 
of duality for the spin chirality. 

2.2 phase modulation 

We consider the Hamiltonian (2) again by introduc- 
ing a phase modulation for the XY term. Such an ex- 
tra phase factor is often discussed in the context of the 
persistent current associated with the gauge flux for the 
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periodic ring. Here, we assign the modulation only for 
the nearest-neighbor bonds, such that 

ee Ji£^(e^+S- +1 + e -^-S+ +1 ) 

n 

+AS*S* +1 ] , (8) 

where is a phase variable. Then the operator of the 
vector chirality order parameter can be obtained as 



- 1 V" - 1 



3 X L d4> 



(9) 



i=0 



where L is the length of the chain. Writing the ground 
state energy per bond with the modulation (j> as E{4>), 
we can then calculate the chirality order parameter 

1 dE{4>) 



Ji 



(10) 



indicating that <f) can be regarded as the field conjugate 
to the chirality order parameter. We can thus discuss 
E((f) and (j>, similarly to the relation between the free 
energy and the magnetic field for the Ising model; if the 
spontaneous chiral order that breaks the Z 2 symmetry 
occurs, it can be extracted as a singularity of E((j)) in 
the limit (j> — > ±0. 

Here, we want to comment on the implication of 
the phase modulation to the experiment; the imaginary 
part of the phase factor in cq. (8) can be casted as the 
Dzyaloshinskii-Moriya (DM) interaction, 17 ) by renormal- 
izing the coupling J cos <j) — > J and J sin (p — > D z , where 
D z is the z component of the DM vector. In this sense, 
it is possible to control the chiral order through the DM 
interaction in the experimental situation. 

2.3 DMRG 

In order to investigate the transformed Hamiltonian, 
we use DMRG. 12 ) Although numerical variables become 
complex numbers, the Hamiltonian is of course an Her- 
mite matrix, so that we have no essential difficulty in 
DMRG computation. For the bulk properties of the chi- 
ral order, we particularly employ the product wavefunc- 
tion rcnormalization group, 13 ) which is an infinite-size 
algorithm of DMRG assisted by the recursion relation 
for the wavefunction. After sufficient iterations, we cal- 
culate the expectation value of the chirality operator in 
the bulk limit as a solution of the self-consistent equation 
of the DMRG recursion relation; since the Z 2 symmetry 
is directly represented as a sign in eq. (6), one of the 
broken symmetry states can be selected in the iteration 
of DMRG, owing to the perturbation from the numerical 
round-off error, or an initial condition of the DMRG. We 
also perform the finite system size calculation of DMRG 
for a finite but long chain, for which the broken symme- 
try state is eventually achieved owing to the numerical 
perturbation as well. We then discuss the spatial distri- 
bution of the chirality. 



Here we should mention about the eigenvalue of the 
density matrix in DMRG. For the chiral-ordered state 
of the original Hamiltonian (4), the density-matrix spec- 
trum has a trivial two-hold degeneracy, where the super- 
position of the two equivalent state originating from the 
Z 2 symmetry is observed. For the transformed Hamil- 
tonian (6), however, we have found that such two-hold 
degeneracy disappears in the following DMRG computa- 
tion. This implies that one of the broken symmetry states 
can be extracted properly. 

3. zigzag XY chain 

For the zigzag XY chain, the chiral-ordered state 
was analyzed by the bosonization 2 -* and the ground state 
phase diagram for the chiral order was obtained in Rcf. 
p], on the basis of the long-distance dumping of the 
chirality correlation function. The phase boundary is 
a c ~ 0.794; the chiral order with the gapless magnetic 
excitation appears for a < 0.794, while the dimer gapped 
state without the chiral order is located in a > 0.794. 18 ) 

3. 1 chirality order parameter 

Since we use the infinite DMRG, the bulk values of 
the chirality is measured as an expectation value of the 
local chirality operator at the center of the chain. In the 
following, we thus denote the bulk expectation values 
of the chiralities for the nearest-neighboring spins and 
ncxt-nearest-neighboring spins as 



Kl EE (K„, n+ i) and K 2 EE (K n , n+ 2>, 



(11) 



respectively. Figure 2 shows the results for k\ and k 2 . We 
can obtain sufficiently reliable results within the number 
of retained bases m — 600 and the number of iterations 
1200, except in the very vicinity of a c and in the region 
a < 0.24. The value of K\ can be checked to be consistent 
with that obtained from the long-distance behavior of the 
chirality correlation function. Note that much more bases 
is required for a < 0.24, since DMRG is not efficient in 
the nearly decoupling limit of double chains. 

In the figure, we can see that the chirality order pa- 
rameter K\ develops almost linearly with respect to a 
for a < 0.5, where we also plot k = 1.73a for compar- 
ison. The scale of the amplitude of the chiral order is 
basically proportional to the zigzag coupling J\. As a 
increases further, the amplitude of K\ decreases rapidly 
toward a ~ 0.81 after taking its maximum at around 
a ~ 0.72. Detailed analysis yields this offset of the chi- 
ral order as a — 0.81(see the next paragraph), which is 
basically consistent with a c = 0.794 on the basis of the 
correlation function. 3 ) In addition, we can see that the 
chirality for the next nearest sites n 2 appears in the same 
range of a as K\\ the weak coupling analysis suggests the 
conservation of the global spin current 19 ) 

Ji/ci = -2J 2 k 2 . (12) 

We have confirmed this relation in the DMRG result, 
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Fig. 2. Vector chirality for the XY model. Solid and open circles 
indicate k\ and K2- The broken line indicates k = 1.73«, which 
is a guide for the eyes. 
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Fig. 4. (color online) Spatial distribution of chiralities for the 
zigzag XY model of a = 0.6. The system size is 384 sites. The 
solid and open circles indicate (k n ,n+l) an d (itn,n+2) respec- 
tively. 



indicating that the spin current is circulating in the 
unit triangles and thus there is no net spin current. 
An interesting point in Fig. 2 is that K2 has a curva- 
ture for a < 0.5, in contrast to the linear behavior of 
Ki ~ 0.173a. According to eq. (12), this behavior is con- 
sistent with K2 = — 2aJ2Ki ~ —a 2 . 

Next, let us consider the critical behavior of the chi- 
ral order; for this purpose, we discuss the normalized chi- 
rality K\ = K\/a, because the energy scale of the chiral 
order is basically proportional to a. Since the Z2 symme- 
try is broken, the expected universality of the transition 
is of 2D Ising class, for which the exponent of the order 
parameter is given by j3 = 1/8. We thus plot kf against a 
in Fig. 3. Then, we can verify the linear behavior of K1 8 , 
indicating that the transition is of 2D Ising universality, 
as expected. Moreover, a straightforward extrapolation 
yields an estimation of the critical point, a c — 0.81. This 
value is basically consistent with a c ~ 0.794 in Ref.p], 
but a slight discrepancy can be seen. Here, we should 

[x10" s ] 

4, . . . . , . . . . , , 
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Fig. 3. vs. a plot for the zigzag XY model. The broken line 
indicates the fitting result: kf = 1.80 X 10 _6 (0.810 - a) 



however recall that the double transition — the present 
order-disorder transition of the chirality and the dimer- 
spin liquid transition — is expected to occur at a c within 
numerical accuracy.; 3 -* Thus, improvement in the accu- 
racy in the vicinity of a c is quite difficult in the previ- 
ous approach of the correlation function. We thus think 
that the results of the two approaches are consistent with 
each other. In addition, it should be noted that, if there 
is a coupling between these two modes, the critical phe- 
nomenon itself may be modified from the Ising univer- 
sality. 

In Fig. 4, we show the spatial distribution of the 
chirality on a finite-length chain of 384 spins, which is 
calculated by the finite-size DMRG up to m = 250. The 
chirality is suppressed at the edges of the chain, but de- 
velops toward the center. In the flat region around the 
center of the chain (n ~ 192), we have (k ntn +i) — 0.0927 
and (k n .n+2) — 0.0278, which are consistent with the 
bulk values in Fig. 2. In addition, we can see that the 
boundary effect more significantly appears for (/t n ,n+2}j 
which splits into two branches near the edges depending 
on even and odd sites. However, it can be confirmed that 
the conservation of the local spin currents 

J\\K>n-l,n) ~ «^2\Kn-2,n) = Jl (&n-l,n) ~ JiV^n^+i) (13) 

is satisfied at every site in the chain. 

3.2 phase modulation 

We next consider the effect of the phase modulation 
in eq. (8). In particular, we focus on the (/(-dependence of 
the groundstate energy and the vector chirality around 
= 0. Using the infinite system size DMRG, we com- 
pute the gourdstate energy per bond E{(f>) and the chiral- 
ity order parameter ki(4>) in the bulk limit. In practical 
DMRG computation, we use matrices obtained for a pre- 
vious (f> as an initial condition for shifted slightly. Thus, 
we have a kind of "hysteresis" in the DMRG calculation, 



J. Phys. Soc. Jpn. 



Full Paper 



-0.327- 




-0.342- 



0.1 
£" 
-0.1 
0.1 

y o 

-0.1 







: a=0.7 




' — i — i — i — i — i — i — 


— I — I — I — I — I — I — : 


; a=1.0 


: 

o 

o 




o 

o 


o 



o 

o 




o 






-0.005 



0.005 



Fig. 5. (a) E(cj>) for a = 0.7(solid circles) and 1.0 (open circles), 
(b) Ki(0) for a = 0.7 (solid circles) and 1.0 (open circles). The 
vertical broken line indicates d> = 0. 



which enables us to deal with a metastable solution in 
the vicinity of <ft = 0. 

Figure 5(a) shows the groundstate energy per bond 
E{4>) and Fig. 5(b) shows the chirality K\(4>) for a = 0.7 
and 1.0. For a = 0.7, the energy has a cusp singularity 
at <fi = 0, while for a = 1.0 the curve is smooth at <j> = 
0. Accordingly, n x for a. = 0.7 has a discontinuity at 
cf) = 0, while that for a — 1.0 is a continuous function 
of <j>. The overhang branch at around <j> = for a = 0.7 
indicates the metastable states associated with the first- 
order transition with respect to <f>. Note that K\ at <f> = 
is agree with that obtained for the unitary-transformed 
Hamiltonian (6) in the previous subsection. It should also 
be noted that the numerical derivative of E(<j>) in Fig. 
5(a) is consistent with J\K\ in Fig. 5(b). This implies 
that the chiral order is realized for a < a c and that <f> 
plays a role of the Z 2 symmetry breaking external field 
conjugate to the chirality. 



4. zigzag Heisenberg chain in magnetic field 

The magnetic phase diagram of the zigzag Heisen- 
berg chain was obtained in Ref. [ 15 ^]. However, the chiral 
order is not detected by the uniform magnetization mea- 
surement. The bosonization study of the zigzag chain 
reveals that the chiral order can emerge in the anti- 
symmetric sector of double chains, rather than the two- 
component TL liquid. 5 ^ Recently, the correlation func- 
tion analysis presents that the chiral order actually ex- 
ists in certain ranges of magnetization. 6,7 ) Moreover, the 
phase boundary between the two-component TL liquid 
phase and the chiral phase can be suggested in the two- 
component TL liquid region of the previous phase di- 
agram, although the magnetization curve has no clear 
singularity. In this section, we would like to address this 
problem on the basis of the direct calculation of the chi- 
rality order parameter based on eq. (6). 

In order to deal with the chiral order in a magnetic 
field, wc employ the finite system method of DMRG, 
since the infinite system algorithm is not suitable for 
dealing with a magnetization-fixed system. We calculate 
the spatial distribution of the local chirality order pa- 
rameters rather than the chirality correlation function. 

4-1 distribution of chirality 

In Fig. 6(a), we illustrate the spatial distribution 
of the local chirality order parameter as a function of 
the site index n and the magnetization M, for a 192- 
site system of a = 1.0. A line in the upper layer in 
Fig. 6(a) indicates (/c n ,„+i) for each magnetization, and 
the density plot at the bottom represents the projec- 
tion of the amplitudes on the n-M plane. We can sec 
that the chiral order appears in the low-magnetization 
rcgion(0 < M < 0.05) and in the upper-middle magneti- 
zation region (0.24 < M < 0.365) ; 20) although the chiral 
order is suppressed at the edges, it develops toward the 
center of the chain. Note that the amplitude of the chi- 
ral order for 0.24 < M < 0.365 is larger than that in 
< M < 0.05. 

In comparison with the previous phase dia- 
gram 15 ) (see also Fig. 7), the region < M < 0.05 cor- 
responds to the field range between the critical field H c 
and the low field cusp of the magnetization curve, where 
H c is the field that collapses the dimer gap. The region 
0.05 < M < 0.24 corresponds to the even-odd (EO) 
phase around the 1/3 plateau, which can be explained 
by the non-chiral single component TL liquid attributed 
to the bound state of the domain wall excitations. 16 ) In 
the region M > 0.24, the chiral order phase clearly ap- 
pears, although the previous phase diagram would have 
indicated a two-component TL liquid. Above M ~ 0.365, 
however, the chiral order disappears again, in spite of the 
lack of explicit anomaly in the magnetization curve. This 
result is consistent with the correlation function analysis 
in Ref. p], suggesting that there is a new phase bound- 
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Fig. 6. (color online) (a) Distribution of the amplitude of the vec- 
tor spin chirality {re n ,n+i} and (b) distribution of the local mo- 
ment (S%). The horizontal axis indicates the magnetization of 
the system M, and the site index n. The vertical axis represents 
k in (a) and (S 1 ^) in (b). In (a), the line in the upper part indi- 
cates («n,n+l) for each magnetization, and the density plot at 
the bottom represents the projection of the amplitudes on the 
n-M plane. 



the saturated state. 14 ) 

At the boundary between the chiral order phase and 
the two-component TL liquid phase, there is no clear 
anomaly in the magnetization curve. However, the am- 
plitude of the chirality order parameter seems to change 
very steeply from a finite value to zero, which suggests 
that the transition might be of first order. Of course, 
the present calculation is for a finite system size, where 
M takes only discrete values. Thus, further research is 
clearly needed for the determination of the transition be- 
tween the chiral phase and the two-component TL liquid. 
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Fig. 7. (color online) Phase diagram of the zigzag Heisenberg 
chain in a magnetic field. The chiral-ordered phase is added to 
the previous one in Ref.[ 15 ']. P: 1/3 plateau, D: dimer-ordered 
phase of zero magnetization, EO: even-odd phase where the os- 
cillating behavior appears in the magnetization curve of finite 
size, depending on the total S z =even or odd. TL1: single com- 
ponent TL liquid, TL2: two-component TL liquid, and chiral: 
chiral-ordered phase. " x " symbols denote the boundaries be- 
tween the chiral-ordered phases and two-component TL phases. 



ary of the chiral-ordered phase at M ~ 0.365. 

In Fig. 6(b), we show the distribution of the local 
spin moment (S%). An important point is that the finite- 
size oscillation of (S*) is very suppressed in the chiral- 
ordered region, and a small but clear 2Uf oscillation of 
the local spin moment (S^) can be seen, where 2fcp is 
the Fermi wave number depending on the magnetiza- 
tion of the system. Since the bosonic field representa- 
tion of Kj, n +i contains a 2kp oscillating term in higher- 
order contributions, the 2k f oscillation appears in the 
chirality-chirality correlation function. 6 ) The present re- 
sults verifies that the same 2kp oscillation can be seen in 
(k ny7l+ i) and (S£). On the other hand, for M > 3.65, the 
oscillation of (S*) exhibits a complicated behavior; this 
suggests that the two-component TL liquid phase may 
be realized in 0.365 < M < 1/2, which can be explained 
by the double- well shape of the spin wave dispersion from 



4-2 phase diagram 

We have performed DMRG calculations for 192-site 
systems of various a and M, and then summarize the re- 
sults in Fig. 7: the magnetic phase diagram of the zigzag 
Heisenberg chain, by adding the present results for the 
chiral order into the previous magnetic phase diagram 
in Ref. [ 15 ^]. The chiral order phase appears below the 
low-field cusp for J2/J1 > 0.7, and above the EO phase 
for J2/J1 > 0.83. The boundaries of the chiral-ordered 
phases and the two-component TL phases are denoted by 
" x " . The boundary of the upper chiral phase approaches 
J2/J1 ~ 0.83, at which the magnetization cusps merge 
into the 1/3 plateau. As J%j J\ increases, this boundary 
rises rapidly toward the saturation field. Also below the 
1/3 plateau, the boundary of the lower chiral phase and 
the low-field cusp meets the 1/3 plateau at J2/J1 = 0.7. 
As J2/J1 increases, then this boundary approaches the 
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lower critical field corresponding to the dimer gap. An in- 
teresting point is that the low-energy excitation around 
the 1/3 plateau can be described well by the domain 
wall excitations carrying S z = ±1/3 and their bound 
states. 16 ^ Then J2/J1 — 0.7 or 0.83 are located at points 
where the nature of the dominant excitations switches 
between the single domain wall and the domain wall 
bound state. This suggests that the occurrence of the 
chiral order can be understood in the context of the na- 
ture of domain wall particles, which may provide another 
viewpoint on the mechanism of chiral ordering. 

5. summary and discussion 

We have discussed the vector chiral order for zigzag 
spin chains. We have introduced a unitary transforma- 
tion that rotates spins on even sites by angle ±7r/2 
around the z-axis. This transformation makes the spin 
current nature of the phase difference more explicit and 
enables the direct numerical calculation of the chirality. 
We also discuss the phase modulation as a conjugate field 
to the chirality. 

For the S = 1/2 zigzag XY chain, we have actually 
calculated the chirality order parameter in the bulk limit, 
using infinite-size DMRG. The result has demonstrated 
that the chiral order emerges for a < 0.81, which is con- 
sistent with the known results. However, we note that the 
unified analysis not only for the transitions of the chiral 
ordering but also for that of the dimer order is still an 
important issue. For the 5 = 1/2 zigzag Heisenberg chain 
in a magnetic field, we have also calculated the chirality 
profile of 192-site systems. We then find that the chi- 
ral order is certainly realized above or below the cusp 
singularities of the magnetization curve, which quanti- 
tatively agree with the correlation function analysis. 6, 7 ) 
The chirality result also suggests that the transition be- 
tween the chiral-ordered phase and the two-component 
TL liquid occurs, although there is no clear singularity in 
the magnetization curve. We have finally illustrated the 
magnetic phase diagram of the zigzag Heisenberg chain 
in Fig. 7. The transition between the chiral phase and the 
two-component TL liquid seems to be of the first order 
within the 192-site result, accompanying a finite jump of 
the chirality. However, further researches of the critical 
phenomena associated with the chiral order are clearly 
needed. 

In order to discuss the physical background of the 
phase diagram, it may be important to see that the cusps 
and transition lines of the chiral order merge into the 
1/3 plateau at a ~ 0.7 or 0.82. The continuous field 
theory states that the chiral order is induced by the Zi 
symmetry breaking in the antisymmetric sector of two 
bosons originating from the double-chain structure of the 
zigzag lattice. On the other hand, the domain wall excita- 
tion carrying S z = ±1/3 also describes the single chain- 
double chain crossover around the 1/3 plateau; 16 ^ The 
nature of the low-energy excitations changes between the 



single domain wall and the bound state of the domain 
walls at a ~ 0.7 or 0.82. Since the cusps and transition 
lines of the chiral order meet at the same point, the do- 
main wall picture and the chiral order should be closely 
related to each other. Such an analysis of the connection 
between chiral ordering and domain wall excitations may 
be an interesting problem. Then, we should also recall 
that the crossover between the interchain and intrachain 
processes of the excitations in the two-component TL liq- 
uid was discussed for the integrable model, which has an 
explicit chirality term in its Hamiltonian. 21 ' 

Finally, we would like to note that the present 
method can work effectively in the ferromagnetic region 
of the nearest-neighbor coupling J\ < 0, where a similar 
chiral order can be expected. 1 ' 22 > 23 ) Although the zigzag 
chain is a simple model, we think that it still contains 
important intrinsic physics. 
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